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^ | Abstract 

O ■ 

■ A new approach to nonperturbative calculations in quantum electrodynamics is pro- 

posed. The approach is based on a regular iteration scheme for solution of Schwinger- 
| Dyson equations for generating functional of Green functions. The approach allows one 

to take into account the gauge invariance conditions (Ward identities) and to perform 
, the renormalization program. 

t^j- | The iteration scheme can be realized in two versions. The first one (" perturbative 

vacuum") corresponds to chain summation in the diagram language. In this version 
the exact result of two-dimensional Schwinger model is reproduced at the first step 
of calculations, but in four-dimensional theory the non-physical singularity (Landau 
pole) arises which leads to the triviality of the renormalized theory. The second ver- 
sion ("nonperturbative vacuum") corresponds to ladder summation and permits one 
to make nonperturbative calculations of physical quantities in spite of the triviality 
i ' problem. 

For chiral-symmetrical leading approximation two terms of the expansion of the 
first-step vertex function over photon momentum are calculated. The formula = 
a/(2ir — a) for anomalous magnetic moment is obtained (a is the fine structure con- 
stant). 

For linearized equation of leading approximation a problem of dynamical chiral 
symmetry breaking is considered, the calculations are performed for renormalized the- 
ory in Minkowski space. In the strong coupling region a > ir/3 the results correspond 
to earlier investigations performed in Euclidean theory with cutoff: solutions arise with 
breakdown of chiral symmetry. For the renormalized theory a solution with breakdown 
of chiral symmetry is also possible in the weak coupling region a < vr/3, but with a 
subsidiary condition on the value of a which follows from the gauge invariance. 
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Introduction 



The problem of nonperturbative calculations in quantum electrodynamics (QED) arose prac- 
tically simultaneously with the principal solution of the problem of perturbative calculations 
which based on renormalized coupling constant perturbation theory. It is necessary to rec- 
ognize, however, that the progress in the nonperturbative calculations during last decades 
is not too large. Quantitative description of nonperturbative effects either is based on non- 
relativistic foundations ( an example is the bound state description based on non-relativistic 
Coulomb problem) or is rather open to injury for a criticism. Besides, the problem of inner 
inconsistency of QED exists (see M). This problem can be formulated as a deep-rooted the- 
sis on triviality of QED in the nonperturbative region (see, for example 0, || and references 
therein) . 

The triviality means that the only non-contradictory value of the renormalized coupling 
is zero. Absence of asymptotic freedom in QED and unsuccessful looking for another type 
of self-consistent ultraviolet behavior are strong arguments in favor of the triviality^. An ex- 
tremal expression of this point of view is a statement that QED can be treated exceptionally 
as the renormalized coupling constant perturbation theory, and, consequently, any nonper- 
turbative calculations are excluded from the consideration. At a high account such a situation 
is not desperate since at very high energies the QED becomes a part of a grand unification 
theory which is based on non-Abelian gauge theory with the self-consistent asymptotically 
free ultraviolet behavior. However, it sounds rather unnatural that for a consistent relativis- 
tic calculation of positronium energy levels, for example, one should exploit a grand unifica- 
tion theory, while it is not need for calculations of annihilation cross-sections. It is difficult 
to find any physical reasoning for such a principal difference of these problems. Therefore 
a possibility to perform approximate nonperturbative calculations of physical quantities in 
the framework of QED itself without contradictions with the triviality of the exact theory ( 
as it made in the perturbation theory) seems to be a necessary component of the theory. 

At first sight the main problem of nonperturbative calculations in QED is an absence of a 
universal small parameter besides the fine structure constant. Due to this reason any partial 
summation of perturbative series seems to be an arbitrary procedure which can be apologized 
only by physical meaning of results. At the same time the absence of a small parameter itself 
is not an obstacle for using some approximation. A motivated approximation can give quite 
satisfactory results even without an explicit small parameter: examples are applications of 
variational methods or mean field type approximations in different areas of physical theory. 

A general problem of various nonperturbative approximations in QED is consistent taking 
into account of requirements of the gauge invariance and renormalizability. It is clear that a 
necessary condition for such taking into account is the existence of a regular iteration scheme 
which in principle permits one to do an arbitrary large amount of steps toward an exact 
solution of the problem. In present work a such type scheme is proposed. The scheme is based 
on an approximation of Schwinger-Dyson equations (SDEs) for the generating functional of 
Green functions of QED by an exactly soluble equation. Its solution generates a linear 
iteration scheme each step of which is described by a closed system of integro-differential 
equations. The requirements of gauge invariance (in the form of Ward identities) are easily 
taken into account at each step of iterations. The renormalization of equations of each 
step is also not a principal problem. Note that equations for Green functions at leading 
approximation and at the first step of iterations look as familiar ones. Similar equations 

In QED with dynamical chiral symmetry breaking the situation can be different (see @,§). 
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were written and investigated earlier in other contexts. A new thing is its appearance in the 
structure of the regular iteration scheme, and it is this circumstance that allows one to give 
them a successive quantum-field-theoretical interpretation. 

Using a bilocal fermion source gives a possibility to formulate the iteration scheme in two 
versions. First of them on the language of Feynman diagrams of perturbation theory is analog 
of the summation of chain diagrams with fermion loop. This version is named " calculations 
over perturbative vacuum" since a unique connected Green function of the leading (vacuum) 
approximation is the free electron propagator. The calculations over perturbative vacuum 
for two-dimensional electrodynamics give the exact result of Schwinger model for a photon 
propagator as early as at the first step, but for physical four-dimensional case they lead to 
the appearance a non-physical Landau pole in the photon propagator and, as a consequence, 
to the triviality of the renormalized theory. Thus a practical meaning of the calculations 
over perturbative vacuum for four-dimensional QED becomes equal zero, though they grasp 
the main nonperturbative effect — the triviality of the full theory. 

The second version of the iteration scheme can be compared on the diagram language 
with a ladder summation. This version is named " calculations over nonperturbative vacuum" 
since the electron propagator of the leading vacuum approximation is a solution of a non- 
trivial nonlinear equation. For this version of the iteration scheme the nonperturbative 
calculations become possible without contradiction to the triviality of full theory. The basic 
part of the work is devoted to the investigation of this version. Apart from the formulation 
of general states and principles of renormalization of Green functions, the article contains 
the calculation of the first-step vertex function (in the case of chiral-symmetric solution of 
the leading approximation equation). Also, for a linearized version of the theory the problem 
of dynamical chiral symmetry breaking is investigated. 

A structure of the paper is the following: in the first Section the necessary notations and 
definitions are given; SDEs for the generating functional of Green functions are introduced 
in the formalism of bilocal fermion source, and generating relations for Ward identities 
are considered. In the second Section a general construction of the iteration scheme for 
solution of SDEs is given; also the leading approximation and the first step of iteration 
over perturbative vacuum are considered. In Section 3 the renormalization of the first-step 
equations is made; it is demonstrated that for the four- dimensional theory the calculations 
over perturbative vacuum lead to the trivial theory already at the first step. In Section 4 
a scheme of calculations over nonperturbative vacuum is formulated. In Section 5 two first 
terms of expansion of the vertex function in photon momentum are calculated for chiral- 
symmetric vacuum. These calculations allow one to obtain a simple formula for anomalous 
magnetic moment: j'2 = a/(2ir — a), where a is the fine structure constant. Section 6 is 
devoted to investigation of dynamical chiral symmetry breaking in QED. In Conclusion a 
discussion of the results is given. 

1 Schwinger-Dyson equations and Ward identities 

We shall consider a theory of a spinor field ip(x) (electron) interacting with an Abelian 

gauge field A^x) (photon) in n-dimensional Minkowski space with the metrics 

xl — x\ — ■ ■ ■ — x\_ x . (For notation simplicity we write all vector indices as low ones) The 
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Lagrangian with a gauge fixing term has the form 

C = -\f^ v - ^(<W) 2 + ~ m + eAty. (1) 

Here = d^A v — d u A^, A = A^j^, ijj = V'*7o, m is an electron mass, e is a charge 
(coupling constant), di is a gauge parameter, 7^ are Dirac matrices. 

A generating functional of Green functions (vacuum expectations values for T-product 
of fields) can be represented as a functional integral 

G(J, V ) = J D(^,A)expi{ J dx(C + J,(x)A^x)) - J dxdy^{y)^ a {y,x)r{x)}. (2) 

Here J M (x) is a source of the gauge field, and r)@ a (y, x) is a bilocal source of the spinor 
field (a and (5 are spinor indices). Normalization constant omitted. 

Functional derivatives of G with respect to sources are vacuum expectation values 

;K/ • < I A^x) I >, s _£f =t<0 I T{r(xW(y)} | > . (3) 

\V ■> x) 



SJJx) ' ^ ' ' ' Srf 



A heuristic derivation of SDEs for the generating functional G is based on relations (see, for 
example, §, 0) 

D{ip,^,A) — expij y + J M (x)A M (x)) dxdy${y)r}(y,x)i(>(x)}, (4) 

= y D(V>, ^, A) ^ 4>(y) expi{ ^ da;(£ + J^x)A^(x)) - J dxdy^(y)r](y, x)ip(x)]. (5) 

Taking derivatives of eqs. (Q)-(^J) and taking into account eq.([|) we get SDEs for the gener- 
ating functional of Green functions of QED 

(9»ud 2 - d,d u + 1^)--^ + zetr^-^— } + J,(x)G = 0, (6) 
di ioJ u (x) L orj(x, X) J 

ev /-A , SG e 5 2 G fit/ is $G . . 

6{x - y)G + (id - m) — + -7^ - - / dxr/ x,a; — = 0. 7 

oij{y,x) 1 dJ^(x)dr}(y,x) J or](y,x') 

(Here and everywhere below denote a differentiation with respect to variable x, while a 
differentiation with respect to another variable will be denoted by indication of this variable 
as an upper index. For example, the differentiation in variable y will be denoted as d^.) 

We consider in this and following Sections the unrenormalized theory, therefore for all 
singular divergent quantities some regularization is implied. A renormalization of SDEs is 
considered below (Section 3). 

Gauge invariance imposes some limitations on solutions of SDEs (^) and (0) known as 
Ward identities. Acting by on eq.@, we get a relation 

^d%--^- + etr Ud-^-\ + dMx)G = 0, (8) 
d\ ibJyXx) L dr](x, x)> 

which we shall name the first generating relation. One more relation follows from eq. ([?D and 
its conjugated SDE 

5{x- y)G + 5G A -id y - m) + - i \f r7 M - / dx' 5G r](x', y) = 0. (9) 



4 



Eq.(^|) is a consequence of the relation 

= J D{^,iP,A)j—rTp{x)eicpi{ J dxiC + J^A^x))- J dxdy${y)r){y,x)ip{x)), (10) 

which is conjugated to (^). Subtract eq.(^) from eq. (]?]), take a trace over spinor indices and 
put y = x. As a result of these simple manipulations we get the relation 

tr{id -, A = / dx'tr{7](x,x') j— - — ry(x',x)}, (11) 

L dr]{x,x) J J 1 dr][x,x') or){x',x) > 

which we shall name the second generating relation. Relation flTT|) is fulfilled really for 
numerous class of interactions local with respect to fermions (see f8|). For QED we can 
combine both relation into one, substituting tr \id s ^u x \ } f rom (0) (i)- As a result we 

get 

i d2d »jW) =dMx)G+e J^ (12) 

which is the generating relation of Ward identities for QED. 

Differentiating eq.(|Ti2|) with respect to J\ and switching off the sources we obtain the 
known relation 

d 2 d ll D llX (x-y) = d l d x 5(x-y), (13) 

where 

6 2 G 



D^\(x - y) 



(14) 

J=7]=0 



5J x (y)5J^(x 

is the photon propagator. From relation ( |I~3"D it follows that the longitudinal part of the full 
photon propagator (in momentum space) is as follows: 

D lT(k) = -d^. (15) 

Differentiating eq.([l|) with respect to r], we obtain after switching off the sources another 
known relation 

^rd 2 d,F^x- x\ y') = e[5(x - y') - 5{x - x')]S{x' - y'). (16) 



Here 

F (ar x' y') = - — 

M ' ' 5J^x)5r)(y',x') 

is the three-point function, and 

5G 



J=T]=0 



(17) 



S(x -y) = 



5r)(y,x) 

is the full electron propagator. 



J=T]=0 



Relation (pi) looks much familiar if one go to the amputated three-pointer (vertex func- 



tion) T M which is defined as 



r M (x | x 1 , y') = I dxidx' x dy' x S 1 (x' - x\)F u (xi] x[, y[)S 1 (y[ - y')D v ^ {x x - x). (19) 

Exploiting this definition and eq.(^), we obtain the Ward identity for vertex function 

id^T^x | x',y') = e[5(x - x') - S(x - y'^S" 1 ^' - y'). (20) 
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2 Iteration scheme 



At e = SDEs @ and (0) have a solution 

Gfree = exp { J_ jy^ * J y + Tr log(l + £ c * 7?) }, 

where 

^ = id^ud 2 - d^d v + jd^d u ]-\ S c = (m- id)- 1 

are free field propagators, and the sign * denotes a multiplication in operator sense. (Also 
in operator sense we shall treat Tr, as distinct from tr, which corresponds to the trace over 
spinor indices.) The functional G^ ree is a generating functional of Green functions for free 
fields and is a base for an iteration scheme of perturbation theory in coupling constant e. 

For solution of SDEs (|6]) and (^) we shall use another iteration scheme proposed in works 
|| , |L0] . A general idea of this scheme is an approximation of functional-differential SDEs (|6]) 
and (0) by equations with "constant", i.e. independent of the sources J M and r], coefficient. 
Thus we approximate functional-differential SDEs near the point J M = 0, rj = 0. Since the 
objects of calculations are Green functions, i.e., derivatives of G in zero, such approximation 
seems to be quite natural. A circumstance of no small importance is a simplicity of the 
leading approximation of this scheme. Equations for all subsequent approximations also 
can be easy written. Green functions of every order of this scheme are defined as solutions 
of a closed system of equations. Technically this scheme is not much more complicated in 
comparison with the coupling constant perturbation theory, but in contrast to the last one, it 
maintains an information which is inaccessible in any finite order of the perturbation theory, 
i.e. it is a nonperturbative method. 

In correspondence with aforesaid, we choose as leading approximation equations a system 
of functional-differential equations 

(g»»d 2 - d,d u + d,d u )-—— + ietr{^— -} = 0, (21) 

di i bJ u (x) L dr](x, x) > 

Ar'(O) p A2/^(0) 

5(x - y)G<°> + (id - m)—- + ^ , J( ., ( S = 0- (22) 

dr}(y,x) % dJ f ,{x)dr](y,x) 

A solution of leading approximation equations ((2l"l)- (|22|) is the functional 

G {0) =exp{iV M * J^ + TrS^*!?}, (23) 
where V^(x) and S^(x, y) satisfy equations 

(9^d 2 - d,A + ^d u )V v (x) + tetr{^S (0 \x,x)\ = (24) 

5(x-y) + (td-m)S m {x,y) + eV(x)S {0) (x,y) = 0. (25) 

By analogy with the theory of ordinary differential equation we shall call these equations as 
characteristic ones. 

In correspondence with the choice of the leading approximation ith term of the iteration 
expansion of the generating functional 

G = G (0) + + ■ • • + G (i) + ■■■ (26) 
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is a solution of iteration scheme equations 

(9,ud 2 - d,d v + d^)-^- + zetr{ 7 ,-^-} = - J,{x)G^ l \ (27) 

5(x - y)GW + (id - m)-—- + - 7 = / dx' V (x,x')-——. (28) 

drj(y,x) % dJ fl (x)dr](y,x) J orj(y,x') 

A solution of eqs.(^)-(^) is looked for in the form 

G (i) = P (i) G (o). (29) 

On taking into account characteristic equations (|24] ) -(|2~5| ) a system of equations for pW 
assumes the form 

(9,ud 2 - d,d u + d,d u )--—- + ietr{^—, -} = -J,{x)P^ l \ (30) 

, 5P^ e 5 2 P^ ~, . 5P« e o(0) , . <5P« 
(«<9 - m) -=-, r + 77^777^7 r + eV(a:) -=-. r + -7^ ( °> (x, y) - 



5r)(y,x) i 5J^(x)5r](y, x) 5r](y,x) i ' SJ^x) 

SD(!-1) 

dx' V (x,x'){— +S^(x',y)P^ i n. (31) 

Since P^ = 1, it is evident that for any i the functional PW is a polynomial in functional 
variables J and 77. This circumstance is very important since it means the system of equations 
for coefficient functions of this functional to be closed in any order of the iteration scheme. 

This iteration scheme has no explicit small parameter. In some sense, the sources J and 
T) play the role of such a parameter. Expansion (p6|) of the generating functional should 
be treated as an approximation of G(J, rj) near the point J M = 0, rj = 0. In essence, 
instead of the question about the small parameter, one should put a question about the 
convergence of the iteration series. Studying the convergence properties of the expansion 
for the quantum-field-theoretical case is obviously a very complicated problem, and we have 
no any rigorous statement concerning this question. Nevertheless, some qualitative reasons 
can be given in favor of the following supposition: the convergence of such a type expansion 
in any case is not worse than the convergence of the coupling constant perturbation series. 
The reason is following: being the simplest from the practical point of view the coupling 
constant perturbation series is the worst in the mathematical sense. The matter is that a 
small parameter (the coupling constant) is a multiplier at the highest derivative terms of 
the functional-differential SDEs for the generating functional. This means that the equation 
is a singularly perturbed one and, consequently, the coupling constant perturbation series 
is an asymptotic series at best. From the other hand, the proposed expansion is regular in 
the mathematical sense and, consequently, possesses the better convergence properties. (See 
also JlO) for more details.) In addition, note that for a model "zero-dimensional" theory at 
n — 0, for which the functional integral becomes the ordinary ones, and SDEs — ordinary 
differential equations, such an expansion possesses very good convergence properties Jl^ 



Generating relation (0) of Ward identities in the framework of given iteration scheme 
takes the form 



U 2 d u — = d u J v (x)G^ + I dx'tv{r,(x,x')— t --——rjfrx)}, (32) 

di oJ v (x) J L orj(x,x') or)(x',x) > 
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Its corollaries like eqs. Q14]), (pig) and ( |20| ) should be modified correspondingly. These Ward 
identities are very useful tools to control the gauge invariance. As we can see below, the gauge 
invariance requirements, which are expressed in the form of Ward identities, impose some 
strict limitations in given iteration scheme (see, for example, discussions after eq. ( Ptj| ) and 
in Section 5). First step calculations, which are considered in present article, demonstrate 
their compatibility with the gauge invariance (though some question arises for the solutions 
with chiral symmetry breaking, see Sections 6 and Conclusion). Surely, at the moment we 
can say nothing on this problem for higher steps of iterations. 

Consider in more detail the leading approximation described by eqs. and fl25|) . 

It is evident that for a conservation of Poincare-invariance of the theory one should assume 

= 0. Then from eq.(p4|) it follows that 



tr 



{ 7 ^(x,x)} = 0. (33) 



By definition (see (fL8f)) is the electron propagator at the leading approximation, therefore 
it follows from above mentioned Poincare-invariance that S^°'(x,y) = S^°'(x — y), and 

SW(x,x) = SV>(Q)= fj^S^p). 



(27T) 

As was stated above, a regularization is always supposed for all such like expressions. 
As these conditions being fulfilled the solution of eq. fl25f ) is the free propagator 



S<® = S c =(m-idY l . 
For this solution condition ( |3"3"| ) is equivalent to the requirement imposed on a regularization 

J m z — jr 

It is very difficult to imagine an invariant regularization for which this condition would not 
be fulfilled. 

So, the generating functional of leading approximation is 

G (0) =exp{TrS*r/}. (34) 
As follows from eq.(|34]), the unique connected Green function of the leading approxi- 



mation is the electron propagator. Other connected Green functions appear at following 
iteration steps. It is also follows from eq. (|34"D that the generating functional at leading ap- 
proximation does not possess the complete Fermi-symmetry. Really, as follows from the 
definition of generating functional, the Fermi-symmetry implies on full generating functional 
the condition 

S2 ° = 52G (35 ) 

5r)' 3a (y,x)S7]^ a '(y\x > ) Sr)( 3 ' a (y',x)5ri' 3a '(y,x') 

Evidently condition (|35D does not fulfilled for defined by eq . (|34]) . The violation of this 
condition leads particularly to the violation of connected structure of the leading approxi- 
mation two-particle (four-point) electron function which is 



{x, y; x , y ) = 



8rjf }a (y,x)5riP' a '(y',x') 



= S%x-y)S c (x'-y') (36) 

J=r/=0 
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- a term — S c (x — y')S c (x' — y) is missed. 

Such a situation is rather typical for nonperturbative calculational schemes with bilocal 
source (for example, for 1/iV-expansion in the bilocal source formalism), but discrepancy 
of such type are not an obstacle for using these iteration schemes. Really, condition fl55p 
should be satisfied by the full generating functional G which is an exact solution of SDEs 
(Q) and (^). It is clear that an approximate solution may do not possess all properties of 
an exact one. In given case we have just the same situation. Properties of connectivity 
and Fermi-symmetry of higher Green functions, which are not fulfilled at first steps of the 
iteration scheme, restore at subsequent steps. For example, the structure of disconnected 
part of the two-electron function is restored as early as at the first step of the iteration scheme 
(see below). At subsequent steps the correct connected structure of many-electron functions 
and other corollaries of Fermi-symmetry are reconstructed. Such stepwise reconstruction of 
exact solution properties is very natural for the given iteration scheme as it is based on an 
approximation of the generating functional G(J,rj) in vicinity of zero. The Green functions 
are coefficients of the generating functional expansion in the vicinity of zero, therefore only 
the lowest functions are well-described at first steps of the approximation - at the leading 
approximation the electron propagator only. Higher many-particle functions come into the 
play later, at following steps, and relation (|35| ) is fulfilled more and more exactly when we 
go toward exact solution. 

A solution of the first step equation is looked for in the form 

G (D = p(i) G (o) j 

where polynomial is a solution of eqs. ([30"D and ( j31~l) at i — 1 and = 0. This solution 
is 

P« = l - V * Si x) * V + S {1) * V + \ % J»* D$ * J, + J M * Fjp * V + a« * J M . (37) 

In correspondence with above definitions, S-j (x, y; x', y') is a two-electron function, 

(x — y) is a correction to the electron propagator, D^J (x — y) is a photon propagator, 

F\ (x? x , y ) is a three-point function, A\P is a photon-field vacuum expectation value. 

Upper index denotes the first-step iteration scheme quantities. 

Equations for first-step functions follow from eqs.(|30"D,(|31"D at i = 1, = and they 

have the form 

(g^d 2 - d^d u + —d^D^lix -y)- ietr {^^^(y \ x, x)} = g flX S(x - y), (38) 

(g^d 2 - d,d u + jdMFPix | x',y') = etr { 7At 4 1) (x,x;x', 2 /')}, (39) 
(9^d 2 - d^d,, + ^rdrd v )Ap (x) + le tr { ltl S {1) (x, x) } = 0, (40) 

(id - m)sP(x, y- x', y') - ie 1 ,S c (x - y)F^(x \ x', y') = S(x - y')S c (x' - y), (41) 
(id - m)F^(z | x,y) = e^S c (x - y)D$(x - z), (42) 
(id - m)SW(x -y)- ie ltl F^(x \x,y) + e lfl S c (x - y)A^(x) = 0. (43) 
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The system of eqs.(|38|)-([43D at first sight seems to be overfilled: six equations for five func- 
tions. Really one can prove that eq.(|59"D is a corollary of other equations. Further, if one 
imposes on the same condition as for (see fl33|)), then the existence of trivial solution 
A\P = for Aft follows from eq. PD|) . We restrict ourselves to this solution for Aft ignoring 
all others as violating Poincare-invariance of the theory. 
From eq . (f42|) we obtain 



F { ft{z \x,y) = -eJ dx'S c (x - x') lfl S c (x' -y)Dft(x' - z). (44) 
Substituting (01) into (BSJ) , we obtain after simple transformations 

PS)) _1 = W _1 + n^ (45) 

where 

n^(x) = ieHi{ ltl S c {x) lv S c {-x)} (46) 

is a free-electron loop. It follows from Ward identities that IT Ml/ is transversal in the momen- 
tum space: 

n^(fc) = n(fc 2 )vr^, (47) 

where 

is transverse projector. On taking into account eq. (|47|) we finally obtain (in the momentum 
space) 



-k 2 + n(£ 2 ) d '(k 

Then we obtain from eq. ([11]) a two-electron function 



Sft (x, y- x\ y') = -S%x - y')S c (x' - y) + 
+ie 2 [ d Xl dx 2 S c (x - x 1 )^S c (x 1 - y)Dft{x x - x 2 )S c (x' - x 2 ) lu S c (x 2 - y'), (49) 



and from eq. (fEf)- a correction to the electron propagator 

SW( X -y) = ie 2 [ d Xl dx 2 S c (x - x 1 )^S c (x 1 - x 2 )Dft( Xl - x 2 ) lu S c (x 2 - y). (50) 



Thus we have obtained the expressions for all first-step functions. Note, that a discon- 
nected part of function Sft (see ([49D) is the "missing" disconnected part of the leading- 
approximation two-electron function (PB"|). Hence, as was mentioned earlier, the correct 
structure of the disconnected part of the two-electron function is reconstructed at the first 
step of the calculations. At subsequent steps a correct crossing-symmetrical structure of the 
connected part is also reconstructed. 

Before renormalization the formulae obtained are formal expressions due to ultraviolet 
divergences. The question arises: is there an explicit regularization scheme, which is consis- 
tent with the considered iterative expansion? Though in general case we have no positive 
answer for this question for any finite step of iteration, at the same time we have no reasons 
to suppose the usual well-known regularization schemes do not work in our case. In particu- 



lar, the ultraviolet divergence of the free-electron loop (46) can be treated with dimensional 
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regularization. Note, that at n = 2 (two-dimensional electrodynamics) this loop converges 
in dimensional regularization. At m = (Schwinger model) 



U(k 2 



e 

71 



and the photon propagator (in transverse gauge) is 

Dm = 1 



- k 2 



71 



/IV i 



(51) 



(52) 



which coincides with the exact result 11 



3 Renormalization 

So far we have considered unrenormalized theory. To give a physical sense to quantities 
to be calculated in a theory with ultraviolet divergences it is necessary to carry out the 
renormalization procedure. 

The Lagrangian of QED for the renormalized fields has the form 

C r = —F^F^ - ^-{d^f + Z$(id -m r + e r A)i, - 6m$1>. (53) 

Here ip and are the renormalized fields, Z\ and Z 3 are renormalization constants of spinor 
and gauge fields correspondingly, m r and e r are renormalized mass and charge, 5m is the 
counterterm of electron-mass renormalization. We have taken into account the fact, that 
due to the gauge invariance and Ward identities the renormalization constants of the spinor 
field and of the interaction are equal: Z 1 — Z 2 , and the longitudinal part of the gauge field 
is not renormalized. 

SDEs for a generating functional of renormalized Green functions are 

{W - B M + j3^l~ + ie rZl tr + J,(x)G = 0, (54) 



t-/ rv /-A . 5G 5G e r Zi 5 2 G 
o[x — y)G + Z\{io — m r )-r— — om-r— H : — 7 M 



Srj(y,x) Sr/(y,x) i 5J^(x)5r](y,x) 

dx'riix.x')-^, -. (55) 

5r](y,x') 

Generating relation (|I2|) of Ward identities for the renormalized generating functional 



preserves its form under substitution e — > e r . Correspondingly, its corollaries (0), (0) and 
(|2"0f) are preserved with the following distinction: all quantities are renormalized ones. 

Being applied to our iteration scheme, the renormalization procedure means: the generat- 
ing functional G as well as all renormalization constants and counterterms has corresponding 
iteration expansions: 

Zj = Zf + Zf ] + ■■■, 5m = 5m i0 ~> + 5m w + ■ ■ • 
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Thus renormalized equations of leading approximation are 



[zf\g, u d 2 - d,d u ) + LdM^^L + ie T Zf ] tr Kt^t} = 0, (56) 

8(x - y)G® + Z[°\id - m r ) — - 5m® — + — ^ 7 = 0. (57) 

dr](y,x) Or](y,x) % dJ^{x)br](y,x) 



A solution of these equations is the same leading approximation functional ([34]) but now S° 
is the renormalized free propagator 

S c = (m r - id)- 1 . 

At that 

Z® = 1, 5m® = 0, 

and the same condition ( |33| ) is imposed for a regularization. Note that the value of constant 
Z® remains undefined in the framework of leading approximation. This constant is fixed 
at the first step of iteration scheme when the photon propagator comes into play. 

The first-step functional is G^ = P®G®, and taking into account above definitions 
and results of the leading approximations we obtain for pi 1 ) the equations 

[Zf\g»„d 2 - dM + -rd,d v }-—— + ze r tr { 7/1 — -} = -J,(x), (58) 

di ibJ v {x) L or](x,x)J 



[id - m r )— + — 7 M + — 7^6 {x - y) 



5r](y,x) i SJ^S^y.x) i SJ^x) 
= J dx x r]{x, x 1 )S c (x 1 -y) + z[ 1] 5(x - y) + 5m®S c (x - y). (59) 

A solution of these equations has the same form (B7|) with a distinction: the coefficient func- 
tions of polynomial ([37]), which define first-step Green functions, are now the renormalized 
quantities. We do not write out equations for the coefficient functions which coincides with 
above unrenormalized equations (j3~8|)-(|43|) up to evident from ([58]) and (|59|) modifications. 

For three-pointer we obtain the same expression (|4*4]), and for the photon propagator 
in the momentum space we have 

"»' 1 = ^w"-*l (60) 



where II is the same free-electron loop pq)-(|47f) with a substitution e — > e r m — > m r . 

Expression ( fj9f) for two-electron function S® is not changed (with the same substitution 
again). For the first-step correction to the electron propagator we have 



S®(x -y) = ie 2 r J dx 1 dx 2 S c (x - £1)7^(21 - x 2 )D®{x 1 - x 2 )j u S c (x 2 - y) 

+Z®S c {x -y)- 5m m f dx x S c (x - x x )S c [x x - y). (61) 
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If photon propagator (|60|) possesses a pole in the point k 2 = /i 2 which corresponds to 
a particle with mass /i, then for U(k 2 ) the following normalization conditions should be 
satisfied 

rV) = 4V, nV) = 4 0) -i- (62) 

In two-dimensional case (n = 2) 

2 = _£^! /- 1 d — x ^ l ~ x }„ . (63) 

v ; vr Jo m r 2 - x(l - x)k 2 - iO v ; 

In particular at m r = : U(k 2 ) = e 2 /iT, and normalization conditions ( |62|) give 



v 3 

1, /i 2 = e 2 /7r, i.e. the result for Schwinger model is unchanged by the renormalization. 

For four-dimensional theory U(k 2 ) is an ultraviolet-divergent quantity. At the dimen- 
sional regularization (n = 4 — 2e) 

n ( fc2 ) = "f^ (2rr)T(e) /' dx 2 * (1 " ^ (64) 
Z7T Jo (m^ — x(l — xjAr — z0) e 

Here T(x) is gamma- function. In this case the normalization conditions give us the following 
4° , = l-^{;+«l)+log2 I r-log^|}, „' 



■ 2 -0. 



Here a r = e 2 /47r, ^(1) is Euler constant, M 2 is a mass parameter of the dimensional regu- 
larization. 

In the Euclidean region k 2 < the renormalized photon propagator in four-dimensional 
theory possesses a non-physical pole at k 2 ~ — m 2 expj^ 1 } - it is Landau pole JT2J. This pole 
arises also at renormalization-group summation 0, and in the framework of 1/iV-expansion 
13]. A presence of this pole is a serious difficulty and, in particular, prevents one from making 



sensible calculations (or needs some refined construction which seems to be superfluous for 
QED |H). In the coupling constant perturbation theory one can exploit a smallness of a r 
to avoid the troubles, but in used nonperturbative scheme the unique non-contradictory 
possibility is the choice a r = 0, and the theory becomes trivial. 

Therefore, as a result of investigation of the first-step equation of our calculational scheme 
the following conclusion is made: for the four-dimensional QED the renormalization of theory 
leads to triviality as early as at the first step of calculations. This result, from one side, 
demonstrates that proposed scheme correctly catches a nonperturbative content of QED, 
but, from the other side, such a version of the scheme is practically useless in contrast to 
the usual perturbation theory. However, the other version of the same calculational scheme 
exists which permits one to make sensible nonperturbative calculations in QED. To state 
this version is a matter of the remainder of the paper. 



4 Nonperturbative vacuum 

The iteration scheme considered above can be modified in such a manner that it becomes 
insensitive to triviality at least for first nontrivial steps of the calculations, i.e. applicable 
to nonperturbative calculations in QED. It means a possibility to calculate physical quan- 
tities without contradiction with the triviality of theory, in a similar way it is made in the 
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perturbation theory in renormalized coupling constant. Of course, one cannot say that the 
situation is fully analogous with coupling constant perturbation theory. In the framework 
of the perturbation theory one can be sure that calculations in any finite order do not lead 
to pathologies. It would be untimely to allege the same statement for the given iteration 
scheme. One cannot exclude at the moment the possibility of a pathology (which originates 
from the triviality of four-dimensional QED) in a higher step of the iterations. In the case of 
such pathology, the use of a grand unification for its elimination seems to be inevitable. A 
full answer on these principal question can be given only with detail investigations of higher 
steps of iterations. Nevertheless, as we shall see, an analysis of first steps of given iteration 
scheme demonstrates that it includes in a natural way such well-known nonperturbative 
treatments of QED as ladder and rainbow approximations. 

To pass this modification of the iteration scheme (which we shall name calculations over 
nonperturbative vacuum) let us resolve SDE @ with regard to the first derivative of the 
generating functional with respect J M : 



1 5G 



i 8J^(x) 



dx\D c (x 



Xx 



){Ju(x 1 )G 



+ ie tr 7j, 



5G 



5r](xi, x\ 



■}■ 



(65) 



and substitute it into the second SDE (|7|). As a result we obtain the "integrated over A^ 
(in the functional-integral terminology) equation 



8(x — y)G + (id — m) 



5G 



5rj(y,x] 
dxi 



6 



+ — dx x D\ " (x - Xi)7 At — — — -tr^y 
p or]{y,x 



5G 



i 



5G 



5rj(x 1 ,x 1 ) 

+ eD c ^(x - si)J„(sih„ }■ (66) 



Exploiting Fermi-symmetry condition ( |5o] ) let us transform eq.(^) in following manner: 



6(x - y)G + (id 



5G 



m 



+ ie' 



5rj(y,x) 
dxArj(x, xi 



dx\D° (x 



5G 



5G 



5r](y,xi 



5r](xi,x) 5t](y,xi) 



Sr](y,x) 



(67) 



From the point of view of exact solutions equations (|66|) and (|67D are fully equivalent since 
the transition from eq.(^) to eq.(|67D is, in essence, an identical transformation. However, 
it is not the case for the used iteration scheme since, as was said earlier, condition ( j35[ ) is 
fulfilled only approximately at any finite step of the iteration scheme. Therefore, eqs. ( |6B| ) 
and flBTD lead to different expansions. Eq.(^) leads, in essence, to the same calculational 
scheme as above, but eq. (|67D gives us a fruitful scheme of calculation of physical quantities. 

Before the writing equations of the scheme note that the photon propagator is not con- 
tained directly in coefficient functions of polynomials of the iteration scheme for eqs. 
( [66] ) and (|6T| ) since these equations are, as pointed above, the result of "integration" over 
the photon variable A^. For calculation of the photon propagator and concerning quantities 
one should exploit the Dyson formula 



D^x-y) = i 



5 2 G 



5J„(y)5J^(x) 



J=T]=0 



D lu{ x ~y)+ ie J dx 1 D° tip (x-x{)tT'ypF v (y\x 1 ,xi), (68) 



which is a consequence of a single differentiation of eq.(|6 
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In correspondence with the general principle of iteration construction we choose as the 
leading approximation eqs. ( pB| ) and flB7|) with J M = 0, 77 = in their coefficients. These 
equations are 

6(x-y)G^ + (id-m) — - + - dx x D c {x - Xl ) llx — ( rtr 7i ,— - = (69) 

dr]{y,x) % J p dr]{y } x) dr)(xi,Xi) 

and 

S(x - y)G w + (id - m)— + ie 2 / dx x D c (x - x x ) 1{l — - lv — = 0. (70) 

6r](y,x) J * oriix^x) 5rj(y,Xi) 

The functional 

G (0) = exp{Tr5 (0) *r/}, (71) 

is a solution for both these equations, but, whereas the solution of the characteristic equation 
corresponding to eq.(^) is the free propagator = S c with condition (|33"D, a characteristic 
equation corresponding to eq.(|70|) has the form 



[S®]- 1 ^) = (m - id)8(x) - ie 2 D; u (x)j^°\x)^ (72) 

i.e., it is a non-trivial nonlinear equation for S^°\ and so we name this scheme the calculations 
over nonperturbative vacuum as distinct from above scheme of eq.(36) which is based on the 
free solution S c . 

A detailed discussion of eq . (|72|) is given below. At the moment we note that at m = 
(chiral limit) in transverse gauge di = this equation possesses a simple solution 

5(0) = -1/id. (73) 

Really, in coordinate spaceQ at n > 2 

D c (x) = Wl 9 12 , I J t±Eig + (i _ dl )(- - (74) 

In transverse gauge di = the function D^ v (x) possesses an important property ("x- 
transversality" ) 

D%{x) ltl x lv = 0. (75) 

From this property the existence of solution (|73 ) at m = follows immediately. 

Atm 7^ solving of the characteristic equation is rather complicated problem which needs 
application of some approximate or numerical methods. The search of solutions, breaking 
the chiral symmetry, is also of great interest (see below, Section 6). 

An iteration equation for the nonperturbative vacuum in correspondence with eqs. 
and fl70D is 

5(x - y)G {%) + (id - m)— + ie 2 / dx x D c {x - xt)^— -7, 



5rj(y,x) J ^ 5r](x u x) 5rj(y,xi) 

( SG^ 5G [i - X) 1 

dxi\T](x,xi)— + eD c (x - x^J^xx)^— \. (76) 

1 077(2/, xi) M dr](y,x) J 



2 The case n — 2 is discussed below, see Conclusion 
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A solution of the first-step equation is = P^'G^ ', where 



= ~»7 * *v + S {1) *r] + J»* * v . (77) 

With taking into account leading approximation ( |70"|) and characteristic equation ( |72|) we 
obtain for three-point function i* 1 ^ , two-electron function and propagator correction 
the following equations 

Fi 1 \z;x,y) = -e J dx 1 D% l {z - x l )S i ~ \x - x l ) 1 ^\x l - y) + 
+ie 2 f dx 1 dy 1 DUx 1 -y 1 )S^\x-x 1 ) llM Fi 1 \z;x 1 ,y 1 ) lv S ( - \y 1 -y), (78) 



8^ (x, y ; x', y') = -S<°> (x - y')S^ (x' - y) + 
+ie 2 [ dx 1 dy 1 D°(x 1 -y 1 )S(°\x-x 1 )^Si 1 \x 1 ,y 1 ;x\y , ) lv S( \y 1 -y), (79) 



- y) = ie 2 / dx x dy x D c [x x - yi)S {0 \x - xi^sf^xi, y x ; yi, y)^ v + 



+ze 2 / dx 1 ciy 1 ^(x 1 -y 1 )^ )(a;-x 1 )7^«(x 1 -y 1 )7^(y 1 -y). (80) 



The first-step equations d78|) — (p0|) are much more complicated in comparison with the 
first-step equations over perturbative vacuum considered above. In diagram language 
perturbative- vacuum equations (|38l)-(f43D correspond to the summation of "chains", and its 
solutions can be easily written out in a general form (see (H)-(ISn)). Equations (fTST)— (BO) and 
characteristic equation ([F^) in the diagram language correspond to the well-known ladder 
approximation. Such type equations for separate Green functions repeatedly were written 
out and investigated in literature (see |L4|— 1[16|| , and also [17j and references therein) as 



the simplest nonperturbative approximation for exact Dyson equations |T8[ which consti- 
tute an infinite system of engaging equations. In our treatment, in contrast to preceding 
investigations, these equation are not a result of more or less arbitrary truncation of Dyson 
equations, but are a consistent part of the iteration scheme. This circumstance is very impor- 
tant for solving such problems as taking into account requirements of gauge invariance and 
renormalizability, which often became stumbling-blocks at investigation of nonperturbative 
approximations. So, for example, one can often encounter the statement that equation (|72|) 
for electron propagator contradicts Ward identity (|T^) and, consequently is not consistent 
with gauge invariance (see discussion of this problem in [|T7|). However, one should not forget 
that a comparison of an approximate equation (|72|) with the exact Ward identity (|16D is not 
correct, and when an iteration scheme is absent the formulation of this problem, strongly 
speaking, is incorrect itself. In the framework of our iteration scheme this problem is solved 
very simply. As follows from eq. (|32|), Ward identity Fjp in this framework has the form 

^-d 2 d v F^\x;x',y') = e[5(x - y') - 8(x - Ol^V - y')- (81) 
di 



It is easy to prove that for i = 1 relation (|81|) and equation (^) are consistent with the 
first-step equation (|78|) for Ffp, i.e. the requirement of gauge invariance is fulfilled. 
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Turn to the renormalization. A renormalized equation for SG/SJ^ has the form 
1 5G /■ . _ , , „ , r , , s „ „ . 



where 



= - / tfciD^x-X! I ^ 3 ){^(a:i)G ! + Z 1 ie r tr7 1/ — , -} (82) 

J OTjyX\, X\) 



Correspondingly, renormalized SDE fl5"7|) is 



5(x - y)G + Zi(id - m r ) SG - 6m SG + 

or]{y,x) or](y,x) 



5 SG 



+i(e r Z 1 ) 2 / dx x D c {x-xx \ Z^^— „ — 

J dr](xi,x) 5r](y,xi) 

f 6G 6G 

= j dxi ^ xi) H^ + zieDc ^ x - xi 1 z * )uxi)i »w^y (83) 

A renormalized leading-approximation equation is 

S ( x - y )G® + Z[°\id - m f )^ - + 

drj(y,x) orj(y,x) 

d Xl D^(x - Xl | Zi%„— - lu — = 0. (84) 

Equation ( |84j ) has a solution in the form of the same functional ([H]), where S 1 '- ) is a solution 
of renormalized characteristic equation 

[S' (0) ]- 1 ^) = (Zf\m r - id) + 6m®)6(x) - i(e r z[ 0) ) 2 D c ^(x \ Z^j^ {x) lv . (85) 

A renormalized photon propagator is defined by Dyson formula 

D^{x-y) =D°^(x-y \ Z 3 ) + ie r Z 1 J dxxD c w (x - x 1 \ Z^trjpF^x^xx). (86) 
Since F^ = at the leading approximation, then 

n(0) _ jnc /y(0)\ _ 1 /. | ,j 

From the normalization condition for the photon propagator we obtain 



zi 0) = 1, 



therefore, in all above formulae one can replace D c (Z%) with D° . 

In general case the renormalized photon propagator in the momentum space is 

ZV(*) = D%{k | Z z ) + ie r Z 1 D c fip (k \ Z 3 ) J -^-tv { 7 pF„(*;p)}, 

where F v (k;p) is a Fourier-image of the three-pointer: 

F »&*>V) = J (2 t)n (2 tV ^ e-^'^rSip - q + fc)^(fc;p). (87) 
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The transversality follows from Ward identities: 



and, taking into account given definition of function f(k 2 ), we obtain 

D»v(k) = - 1+ JjfK iat - d^-. (88) 



Zero mass normalization condition ("photon normalization ") gives us the following 

Z 3 = l + /(0). (89) 

The constant Z\ and the mass-renormalization counterterm 8m are connected with the 
normalization of the electron propagator. In general case the electron propagator in the 
momentum space is defined by two scalar functions 

n _ i _ Hp 2 ) + a(p 2 )P 



b(p 2 ) — a(p 2 )p b 2 (p 2 ) — a 2 (p 2 )p 2 

If the propagator S(p) possesses a pole at the point p 2 = m 2 ,, then the normalization condi- 
tions in this point are 

b(m 2 r ) = m r a(m 2 ) , (90) 

a{m 2 .) + 2m 2 r a'{m 2 ) - 2m r b'(m 2 ) = 1. (91) 

A chiral limit for the renormalized theory means that the chiral-non-symmetric terms 
disappear from the renormalized Lagrangian 



Z x m T + 5m -> 0. (92) 



For transverse gauge d\ = renormalized characteristic equation ( |85| ) has a chiral-symmetric 
solution at the chiral limit 

S{0) = -i" ( 93 ) 
Z\ >p 

For this solution normalization conditions (BH) and fl9"T|) give us the relations 

Z[ 0) = 1, m r = cW 0) = 0. 

Hence, this solution coincide with unrenormalized solution (|73|). 
Renormalized equations for the first-step coefficient functions are 

F^iz; x, y) = -e r z[ 0) J d Xl D c x ^z - Xl )S^(x - x 1 ) 1 ,S^(x 1 - y) + 
+i{e r zf ] ) 2 J dx 1 dy 1 D^(x 1 - yi )S^(x - x 1 h,Fi 1 \z;x 1 ,y 1 ) lu S^(y 1 - y), (94) 

(x, y;x',y') = -S® (x - y') S® (x' - y) + 
+t(e r z[ 0) ) 2 [ dx 1 dy l Dl u (x 1 -y 1 )S^\x-x 1 )^Si 1 \x 1 ,y 1 ;x',y , )^ \y l -y), (95) 
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S {1 \x - y) = i{e r Zf ] f J dx x dy x D c ^(x x - y x )S^\x - x 1 ) f y lt 3^\x u y 1 \yi,y)nf v + 

+ J dx x S^{x — x\ )\ [Z[ 1] (id - m r ) + 5m^]S {0 \x x -y) + 
+te 2 r Z[ 0) J dytfZjpD^fa - y x ) + Z^D^fa - Vl \ Z^)] llx S^{x x - y x ) lv S^{y x - y)} + 

+i(e r z{ 0) ) 2 f dx x dy x D°{x x - y x )S^{x - x x ) lfl S^{x x - yi)i^\ yi - y).(96) 



,(o) 



As we see, the gauge invariance applies rather strict conditions: constant Z\ is the unique 
renormalization constant for two first-step equations (CT ) and (0). In the equation for 
new renormalization constants appear: there are Z x , Z% (the last one is fixed by condition 
(|89|)) and the first-step mass- renormalization counterterm 5m^. 



5 Vertex 

Let us go from the three-point function F x to the amputated three-pointer — vertex function 
(P~9|) . In correspondence with eq.fl9~4|) the first-step equation for the vertex is 

T x (z; x, y) = -ej x 5(x - y)5(x - z) + 
+ie 2 D c pv (x ~y) J dx x dy x j^,S(x - x x )T x (z; x x , y x )S(y x - y)j v . (97) 

In this Section we shall use simplified notations 

r A = r?\ e = e rZ [°\ s = s®. 

In the momentum space equation ( B7h has the form 

T x (k;p) = -e lx + ie 2 J j^D»u(p ~ p'h,S(p')T x (k; p')S(p' + k) lv . (98) 

Here k is the photon momentum, p is the electron momentum (see definition (R^) of the 
Fourier- image of three-pointer). Equation (]98D is known in literature as Edwards equation 
Note again that in our approach this equation is not a result of arbitrary truncation, but 



is a consequence of the iteration scheme equations. Electron propagator S in this equation is 
a solution of characteristic equation ( jB5f) which defines the leading vacuum approximation. 

We consider a solution of vertex equation fl9~7|) at small A; with above mentioned chiral- 
symmetric propagator ( ]93|) in the transverse gauge = 0. To solve this equation it is 
convenient to introduce a function Q x (k;p) which is defined by the relation 

T x (k;p) = p$ x (k;p)(p + k). 

Expand $a near k = 

<5> x (k;p) = $ A (p) + k p $xp(p) + ■■■ 

Here 

*a(p) = $a(0;p), $ Xp (p) = d k p $ x (k;p)\ k=0 . 

At that 

T Xp (p) = d k p T x (k;p)\ k=0 =p$ Xp (p)p + p$ x (p) lp . (99) 
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Equations for $a(p) and $a p are 

/dp' 

and 



r dv' 

P$\p(p)p = -P^\{P)1 P + ie 2 r J j^D c ^(p-p')j^ Xp (p')j u 



It is more convenient to solve eq. (|100|) in the coordinate space: 

- <9$ A (x)<9 = -e-y x 8(x) + ie 2 T Di v (x)^ p ^ x (x)^ u . 



(100) 



(101) 



(102) 



Here D c (x) is defined by eq.([F3). Due to the property of " x-transversality" (see (|75|)) all 



iterations of eq. (|102|) turn into zero. Indeed, a zero approximation is 



*X(P) = = e<f x \ 

p p p 



and in x-space 

Consequently 
and 



i e -*W2r(n/2) x A x 
2W 2 (x 2 - iO)™/ 2 

$ A = $ °A 



X. 



r A (p) = -e 7 A- (103) 

Let us go to solving eq. (|101|) for $ Ap . (We consider now the four-dimensional case n = 4.) 
This equation is also more convenient for solving in the coordinate space, where it has the 
form 



- d$ Xp (x)d 



e 7a^7p , ,„2nc 



+ <D c (x) 1 ^ Xp (x) lv . 



2tt 2 (x 2 -z0) 2 r "" v 
(We have taken into account the result of solving equation for $ A (p).) Here 

1 



QpLV X p Xy 



Utt 2 (x 2 - iO) 2 x 2 - iO 



(104) 



(105) 



Expanding $ Ap over spinor structures 



$Ap = ®\polo + ®l pa l5lo 



we obtain two equations 



(g ar d 2 - 2d a d T )§ XpcT (x) 



e g Xr x p - g Xp x T + g pr x x 



(g <JT d 2 -2d <J d T )<$>\ pa (x 



2~r? 



2tt 2 
2a r 1 
it x 2 — iO 



X n 



(g a 



{x 2 - zO) 2 

%lC fry 00 *7~ \ 



x 2 — iO' 



2a r 



£\po 



(go 



%JC f-r 00 *r 



(x 2 — iO) 2 it x 2 — iO ar x 2 — iO 



(106) 

)^lpo(x). (107) 
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Here 

e 2 

a r = —. 
4tt 



As for equation for $ A iterations of eq. (|1061 ) for $ Apcr turn into zero, and its solution is 

Xpa[X) = ~4^(x 2 -i0) 2 ' 



Equation ( |107|) for <3> A possesses a solution 



A ^ v ; 47r(27r-a r ) VffT x 2 -zO- 
Going to the momentum space we finally obtain for T\ p (see (ffiD): 

F\p(p) = ie—^- — e Ap<TT 7 5 7 CT Pr (108) 
27T — a r + 2O 

Matrix element of vertex function T\(k;p) at small defines two formfactors fi and f2 

6 

M(g)r A w(p) ~ -e r fiu(q)^ x u(p) - -^hu{q)[^\, k}u{p), (109) 

which connected correspondingly with a physical charge and a magnetic moment (see, for 
example, |H5|). Here u(p) is a solution of Dirac equation (p — m r )u{p) = 0. Comparing with 
calculated vertex we obtain at k = 

e r /i = e = e r z{ 0) . 

For the normalization on the charge e r we put fi = l and, consequently, 

zi 0) = 1. 

This very important fact means that the normalization of the vertex function on the renor- 
malized charge gives in our approach the same value of renormalization constant Z\ as 
wave-function renormalization condition (|9~TD, i.e. our calculational scheme is consistent 
with requirements which imposed by gauge invariance. 

As for the second formfactor f 2 , which defines a correction to the magnetic moment, 
note, that the massless particle, of course, has no magnetic moment. However, the calculated 
first term of vertex function expansion in k^ (see ( |108| )) permits making a nonperturbative 
estimate for anomalous magnetic moment in chiral limit. 

Let us define the problem more exactly. In expansion of f 2 in orders of a r 



7T V 7T ' 



coefficients Cj, starting with the second one, are mass-depended. (So, they have even different 
signs for muon and electron, see, for example, pOfl .) Nevertheless, one can, in principle, to 
consider a peculiar "chiral limit" for these coefficient, i.e. their values at m T — > 0. It is 
this estimate which we are talking about. To obtain this estimate let us make an identical 
transformation of the second term in eq. ( |109|) 



— u{q)h\, k}u{p) = iu(g)[7A, k}pu{p). (110) 

TT\iy P 
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For massive particles it is an identity due to Dirac equation for u(p). However, for r.h.s. 
of eq.( |110[ ) one can use result ( |108[ ) which was obtained for massless particle. Comparing 
coefficients at t\ paT1 we obtain 

(HI) 

Z7T — a r 

The first term of the expansion in a r coincides with the Schwinger correction (see, for exam- 
ple, |2(|). At formal limit a r — > oo the full magnetic moment in correspondence with eq. ( |111|) 
turn to zero. It is interesting to note, that similar screening of chromomagnetic moments of 
quarks was pointed out in some models of nonperturbative quantum chromodynamics [EH 
and also in relativistic quarkonium model |22| . 



6 Dynamical chiral symmetry breaking 

In this Section we consider a problem of dynamical chiral symmetry breaking (DCSB) for 
QED in the framework of proposed calculational scheme over the nonperturbative vacuum. 
Note, that our present consideration does not cover all aspects of this complicated problem, 
in particular, we do not consider a question of connection of DCSB with triviality problem. 
A consideration of this and a number of other questions imposes an investigation of the 
first-step equations for higher Green functions. We limit ourselves to investigation of the 
leading- approximation equation. 

A renormalized equation of leading approximation for the electron propagator has the 
form 

S-\ x ) = (m - Zid)5{x) - ie 2 D c ^(x)^S(x)^. (112) 
Here and everywhere in this Section we denote 

S = S®, Z = Z{°\ m = z[ 0) m r + 5m(°\ e = e r zf\ 

and the free propagator D c (x) is defined by eq.(|73|). 

In general case (if one does not consider parity breaking solutions) S can be represented 

as 

S = id A + B, 

where A and B are scalar functions of unique variable. Correspondingly, for the inverse 
propagator we have 

S- 1 = -ida + b, 

where in the momentum space the functions a and b are connected with the functions A and 
B by relations 

^ = ~L2 2~~2' B — 75 2~~2' (113) 

o — crp z o — a z p z 

Taking into account eq.([Ti[) and properties of Dirac matrices, characteristic equation (|112| ) 
in the coordinate space can be represented as a system of two equations 

e -i7rn/2p/n\ 

ida = Zi05(x) - 2ad lr — ■ idA, (114) 

b = m8(x) + a(l - n - dfi r _^ 2 _ >* ' ■ B. (115) 



7T (X 
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Here 

= i- 

As is evident from eqs. (|114|) -( |iT5|) , two special gauges exist for characteristic equation 

(H) 

1) the transverse Landau gauge di = 0. In this gauge eq. (|114j) becomes trivial and pos- 
sesses a solution a = Z8(x). 

2 ) the gauge d\ — 1—n. In this gauge^ b = mS(x) in case if the product B(x) ■ (x 2 — iO) 1 ^' 1 ^ 2 
is well-defined in a sense of generalized functions, i.e., divergences are absent. Otherwise, 
in contrast with the Landau gauge, we have an uncertainty of the type • oo. Its evaluation 
is performed by a regularization, and we can only say that in this gauge b(p 2 ) = const. 

Below in this Section we shall consider the four-dimensional case. Multiplying eqs.( |114j ) 
and (|115|) by x 2 , we obtain the equations 

2ry 

x 2 ida = diidA, (116) 

7T 

x %= — (3 + dAB. (117) 

7T 

Such a multiplication can be considered as some regularization of potentially singular prod- 
ucts in the r.h.s. of eqs. (|114|) and (|115|) . More, going to p-space and taking into account 
eq.( |113|) , we obtain for a{p 2 ) and b{p 2 ) a system of ordinary differential equations 

d 2 a da a , a , 

t ^ + 3 Tt = ^ d -¥^Ft' (118) 

d 2 b db a , , . b , 

*5? + 2 s = ^ 3 + * ) F^- < 119 > 

Here t = p 2 . 

This very complicated system of nonlinear differential equations has the simplest form in 
the transverse Landau gauge di = 0. Below we shall exploit this gauge. As was mentioned 
earlier, in this gauge a — Z, and system ( |118| )-( pT9|) is reduced to the equation for b 

d 2 b db 3a b ,^ nn . 

t w +2 dt = ^W^zn- (120) 

This equation always has the trivial solution 6 = 0, which corresponds to above mentioned 
chiral-symmetric solution (|93|) . The existence of non-trivial solutions is also possible, and 
their asymptotics is not difficult to find. Really, at t — > oo two variants are possible: 

a) b 2 ~ t; 

b) b 2 <| t \. 

It is not difficult to prove that the third possibility b 2 t \ up to logarithms leads to 
the variant a). 

Consider the variant b). Then in the ultraviolet region t — > oo equation (|120p is reduced 
to Euler equation, and an asymptotic behavior of b(t) is 



b ~ r i/2(i- v /l -3^A) at ar < (i 2 l) 



5 At n = 4 it is Solovev-Yennie gauge EE 
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b ~ r 1/2 \ogt at a r = tt/3 (122) 
and 



J3a r /-K — 1 

6~r 1/2 sin{-^ logt} at a r > tt/3. (123) 

2 

Note, that here a r = |^ is the renormalized fine structure constant. 
The existence of a critical value 

a r = a c = tt/3 (124) 
should be pointed out, at which a change of ultraviolet behavior regime takes place. The 



existence of this critical point was firstly pointed in work | 24fl . In works Q such regime 



transition was connected with DCSB in QED. (See also [I7],[[7[ where this approach is dis- 
cussed in detail and the extensive bibliography is given. For latest developments of this 
approach see [|25|| .) In all cited works a scheme with ultraviolet cutoff in Euclidean space 
was exploited. A main calculational ansatz was a linearization of eq. (|120|) , which consists 
in an approximation F] 

b b 
b 2 — t m 2 — t ' 

in the r.h.s. of eq.( [L20|) . Here m = 6(0). After such a linearization the equation for b(t) be- 
comes a hypergeometrical equation. Boundary conditions for this equation are defined from 
the integral equation which is eq.(|72|) in the momentum space An asymptotic behavior 
of solutions of linearized equation is given by the same formulae ( [[ 2 1| )- (|T23"D , and this fact 



is the main argument in favor of the linearized version (at least in the ultraviolet region). 
These solutions, of course, essentially depend on a cutoff parameter which is included in the 
boundary conditions. Nevertheless, one was succeeded in demonstrating the existence of a 
phase transition at the critical point a c = tt/3. This transition corresponds to DCSB. In the 
chiral limit under the critical point (weak coupling) only the chiral-symmetric solution exists, 
but beyond the critical point (strong coupling), when the ultraviolet asymptotics becomes 
oscillating, a solution with dynamical mass b ^ arises, i.e., the spontaneous breaking of 
chiral symmetry takes a place. 

We consider this problem in the pseudoeuclidean Minkowski space for renormalized equa- 
tion ( |1 12p .pl At that we shall use another approach to the investigation of the propagator 



equation, which, nevertheless, is similar to the approach of above mentioned works and also 



based on a linearization of nonlinear equation ( 112 ). Since we work with renormalized the- 
ory, our results are cutoff- independent. In general they are consistent with the results of 
linearized unrenormalized theory in Euclidean space up to the unique exception: in pre- 
critical region a r < a c DCSB is also possible, but under some condition on value of a r . An 
investigation of this condition requires studying of equation for the three-point function with 
the nonperturbative electron propagator and goes out of the framework of this work. 



We shall solve equation ( |112j) in transverse gauge di = at m ^ by iterations and a 



foundation for iterative solution will be an exact solution at m = 



4 Note, that in the cutoff scheme Z = 1. 

5 An essential point at the formulation of these boundary conditions is exploiting of Euclidean version of 
the theory. 

6 This point seems to be essential since detailed investigations of unrenormalized equation ( |7^ ) in Euclidean 
space at m = (see, for example, |2q] ) demonstrate the existence of non-pole complex singularity for its 
solutions, and, hence, the Euclidean rotation seems to be a problem. 
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At that the linearization consists in rather natural procedure from the point of view of the 
iterative solution: starting from the representation of the inverse propagator 



S- 1 = S X + £ 

we approximate 

S = [V + S]- 1 « S - S * E * S 



Z(p 2 + i0) Z 2 (p 2 + i0) 2 ' 

It is clear that such an approximation can be fully valid only in the ultraviolet region | p 2 
fi 2 , where /1 2 is some scale which can be treated as an infrared cutoff. Since a nonperturbative 
region of QED is the ultraviolet region, the natural supposition is in the following: the 
ultraviolet behavior is crucial for basic nonperturbative effects including DCSB. 

Using the property of x-transversality (see (|75|)) it is easy to prove E = 1 • a, and for 
scalar function a in the coordinate space we obtain the equation 

7T x 2 — iO 

where 



a( x 2 ) = m5{x) + "7 ^ - - ■ <j)(x 2 ), (125) 
a{p 2 



p l + 

The four-dimensional delta-function S(x) in the Minkowski space is a Lorentz- invariant 
distribution (see, for example, |27|), and to solve eq. ( |125| ) it is convenient to use a represen- 
tation of the delta function as a limit of a sequence of functions in x 2 -variable. As such a 
representation we choose the formula 

S 4 (x) = -\\im\(x 2 -i0) x - 2 (126) 

7T 2 A^O 

which can be easy proved with the Fourier-transform of function (x 2 — i0) x ~ 2 (see, for ex- 
ample, p8|). 

A partial solution (a solution of inhomogeneous equation) will be searched as 



a (x 2 ) = Co{x 2 - zO) 



A-2 



at small A. As we shall see below, 1/A plays the role of a regularization parameter. Making 
the corresponding Fourier-transform (see p8||), we obtain 

and, substituting into eq. (|125|) , we finally obtain 

Co = -^-X 2 + 0(X 3 ). 
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Consequently, the solution of the inhomogeneous equation is 



^o(p 2 ) = • (128) 

The full solution of eq.( |125| ) is a sum of o"o and a general solution of homogeneous equa- 
tion. A solution of homogeneous equation is sought for in the form 



a 



(x 2 ) = C(x 2 -i0f- 2 . 



Hence 4>(x 2 ) is defined by the same formula ( |127| ) with the substitution A — > [3. For parameter 
f3 we obtain the equation 

0(1-0) = 2^. (129) 

If one chooses a solution with C = 0, i.e., limits oneself to consideration of partial solution 
( |128| ), then normalization conditions ( |90"D and ([Jl]) give f\ 

5m^ = — rWV> Z — 1 
4vrA 

and, correspondingly, in this case 

S^ 1 = —p + m r . 



In the chiral limit (see (|9~2"D), when m — > 0, we obtain m r = 0, therefore this solution is a 
chiral-symmetric one. 

But if one permits solutions with C ^ 0, the situation is essentially changed. Note at 
once, that the solutions of the homogeneous equation in p-space are singular at p = 0, but, 
due to the essence of the approximation made, all subsequent formulae can be interpreted, 
as was pointed out above, only at | p 2 fi 2 , where fi 2 is an infrared cutoff. 

So, in dependence of value of parameter a r , we have for electron mass function: 

1) At a r < a c (weak coupling) 

b(p 2 ) = -^ + C(p 2 + i 0r^ (130) 



where /3 = |(1 — \Jl — a r /a c ). 
2) At a r = a c = n/3 



b(f) = ^ + ™&L. (131) 



3) At a r > a c (strong coupling) 



3a r (p 2 + i0y/ 2 ' 



4nm\ , Csin{f log^} 



7 The formula for counterterm 8m demonstrates a role of A as a parameter of regularization which removes 
ultraviolet divergence. Therefore, the representation of delta- function by formula (126) can be considered 
as a special analytical reg ular ization 

8 We take a root of eq.(129) which corresponds to the ultraviolet asymptotics of exact solution, see (121). 
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Here 



UJ 



a = Z for any of three cases. 

Below we consider the chiral limit m = Zm r + 5m^ = 0. Solution parameters are fixed 
by normalization conditions (|9l3|)-(|91|) of the electron propagator. If (as in ultraviolet-cutoff 
scheme P|,[17[l, |25|) one takes Z — 1, then the normalization conditions become 



b{m 2 ) = m r , m r b\m 2 ) = 0. 



(133) 



In the weak coupling case (a r < a c ) it follows from the normalization conditions that either 
C = or m r = 0, i.e. in this case DCSB-solutions are absent. In the critical case a r = a c a 
solution with m r ^ is possible which has the form 



b(p 2 ) 



mz 



( P 2 + ioy/ 



/ 1, p 2 + iO\ 

^( 1 + 2 log ^r)- 



(134) 



At a r > a c (strong coupling) a DCSB-solution also exists: 



b(p 2 



mz 



{ P 2 + toy/ 2 



cos 



rCJ, p 2 + i0^ 1 . f W, p 2 + i0> 

l2 log ^r} + ^ sm t2 log ^r} 



(135) 



As we can see, these results correspond to those of the unrenormalized theory with a 
cutoff in the Euclidean space (see [fj]J17||)|25"||). However for the renormalized theory we have 
no any prior reason to set Z = 1. If one refuses this condition, the results in the weak 
coupling region can be essentially different. Really, at Z ^ 1 the normalization conditions 
are 

b(m 2 r ) = m r Z, 2m r b\m 2 ) = Z - 1. (136) 

At a r < a c , in contrast with the case Z = 1, the existence of DCSB-solution becomes 
possible. It has the form 



1 



At a r = a r a solution is 



b{p 2 



l + 2(3 l 



mz 



Zm r (^ 



mz 



p 2 + i0 



(137) 



( P 2 + ioy/ 2 



Z+(Z- - log — 



mi 



:i38i 



and, at last, at a r > a c 



b{p 2 ) 



mz 



( P 2 + ioy/ 2 



Z cos \ - log — \ + -(2Z-l) sin \ - log — \ 

1 2 mi * uo 1 2 mi J 



p 2 + iO - 



(139) 



(a = Z for two last cases, Z is arbitrary.) 

Therefore, with taking into account the renormalization freedom DCSB becomes possible 
in the weak coupling region which corresponds to physical sector of the theory. At that, 
however, an important limitation connected with the gauge invariance arises. The necessity 
for such a limitation follows also from a consideration of switch-off-interaction limit. Indeed, 
it follows from eq. (|137|) that at a r — » b — » m r + 0(a r ), i.e. either m r = 0, or DCSB does 
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not disappear in the switch-off- interaction limit. The last possibility looks quite strange. A 
way out of this situation can be in the following: all the picture in the weak coupling region 
is realized only for some finite values of a r . A condition for such values is a normalization 
of the vertex on physical charge e r : 

T^(k = 0;p 2 = ml \ a r , Z) = -e/y^. (140) 

Due to the gauge invariance the renormalization constant Z is the same as for the electron 
propagator equation, and substituting the value Z = Z(a r ) from eq.( |137| ) we shall obtain 
an equation for a r . Solutions of this equation will define "permitted" values of the charge. 
A realization of this program requires solving an equation for the vertex function with 



propagator ( |137| ). 

At a r > a the situation changes in principle: in this case Z (at the leading approximation) 
is arbitrary, and the normalization on the physical charge fixes Z, but does not lead to 
additional limitations on values of a r . 



Conclusion 

Above results give us a reasoning to consider the calculations over the nonperturbative 
vacuum as an quite adequate scheme for calculations of nonperturbative effects in QED. 
A formulation of system of equations for Green functions at any step of this calculational 
scheme is technically as simple as the coupling constant perturbation theory f\ and does 
not require anything except of the "stupid differentiation". The equations themselves, of 
course, are much more complicated, but, as we can see from the calculation of the vertex 
function for chiral-symmetric solution (Section 5), they possess much more simple solutions in 
comparison with analogous equations for Green functions in nonperturbative approximations 
which based on a more or less arbitrary truncation of Dyson equations. 

As regards the calculations of Section 6, which concern the DCSB problem, the one of 
main questions concerning the linearized equation for electron propagator is a problem of a 
gauge dependence of the results, in particular, of a gauge dependence of the critical constant 
[17| and references therein). In our treatment this equation is not a result of an 



a r see 



arbitrary truncation of Dyson equations but a consistent part of the iteration scheme, and 
this question becomes unavoidable. In this connection we note that the linearization seems 
to be substantiate in Landau gauge only, and for other gauges the nonlinear effects scarcely 
can be taken into account with such simple ansatz. As has been mentioned above, the 
results of Section 6 in strong coupling region coincide in essence to the results of so-called 
rainbow approximation, or "quenched QED with bare vertex" (see [17], ||25|| ). It is no wonder 



since the basic equation for electron propagator is the same, and a difference consists in the 
regularization and renormalization scheme. Generally, quenched QED is an approximation 
of QED with an effect of pair creation being neglected. For the electron propagator this 
approximation permits more general form of the vertex function. For example, the use of 
Curtis- Pennington vertex [39| solves the problem of gauge dependence of critical coupling 



a c , though is not solve all the problems in strong coupling region (see |3(J and references 



9 Note, the coupling constant perturbation theory can be considered as a partial case of the general 
iteration scheme based on principles of Section 4. To obtain a perturbation theory series it is enough 
consider the same iteration scheme with single Grassman fermion sources over perturbative vacuum which 
corresponds to leading approximation G' ' = 1 (see also |9fl-||lC|). 
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therein). Our approach is not an approximation of QED but an iteration scheme, so we have 
not a possibility for such generalization of leading order equation for electron propagator. As 
a general remark, note, that the strong-coupled four- dimensional QED is of pure theoretical 
interest as a model gauge theory for investigation of chiral symmetry breaking phenomenon. 
The results of various approaches (such as worldline variational approach, high orders of 
perturbation theory, quenched approximation etc.) essentially differ in strong coupling region 
(see for review), which indicates our present poor understanding of this problem. 

In closing we shortly discuss a problem of DCSB for quantum chromodynamics (QCD). 
In contrast to QED, where nonperturbative effects (and DCSB among them) are defined by 
ultraviolet region, for asymptotically-free QCD a nonperturbative region is the low-energy 
infrared region. At that a mechanism of DCSB in QCD can be even more simple in compar- 
ison with QED, since effective taking into account of gluon self-action in the infrared region 
inevitably leads to appearance of new dimensional parameters. 

As a simple example of a realization of similar DCSB mechanism by infrared singu- 
larities consider two-dimensional QED. The analogy among the nonperturbative effects of 
two-dimensional QED and those of four-dimensional QCD was repeatedly pointed out. In 
spite of certainly limited nature of such analogy, an investigation of the nonperturbative 
dynamics on example of much more technically simple QED can be useful for understanding 
of QCD in the nonperturbative region. 

At n = 2 the free propagator D c (x) (see (0)) entering equation ( |112| ) in general case 
is infrared-singular, and it is necessary to introduce some infrared cutoff. However, a gauge 
exists for which such a cutoff is unnecessary. This is above mentioned infrared-finite gauge 
di = 1 — n = — 1. As follows from eq. (|74|) , in this gauge 

and the equation for the electron propagator is simply 

S- 1 = (S c )~ 1 + 2aS. (141) 

In the momentum space eq.( |141|) is reduced to a system of algebraic equations, and among 
its solutions a DCSB-solution exists. Surely, in the two-dimensional theory the sponta- 
neous breaking of continuous chiral symmetry is not realized since such a state is unstable 
in correspondence with Mermin-Wagner-Coleman theorem. Nevertheless, the fact itself is 
remarkable from the point of view of above mentioned analogy with four-dimensional QCD. 
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